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Abstract 

The pairing Hamiltonian constitutes an important approximation in many- body systems, it 
is exactly soluble and quantum integrable. On the other hand, the continuum single particle 
level density (CSPLD) contains information about the continuum energy spectrum. The question 
whether one can use the Hamiltonian with constant pairing strength for correlations in the contin- 
uum is still unanswered. In this paper we generalize the Richardson exact solution for the pairing 
Hamiltonian including correlations in the continuum. The resonant and non-resonant continuum 
are included through the CSPLD. The resonant correlations are made explicit by using the Cauchy 
theorem. Low lying states with seniority zero and two are calculated for the even Carbon isotopes. 
We conclude that energy levels can indeed be calculated with constant pairing in the continuum 
using the CSPLD. It is found that the nucleus 24 C is unbound. The real and complex energy 
representation of the continuum is developed and their differences are shown. The trajectory of 
the pair energies in the continuum for the nucleus 28 C is shown. 
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I. INTRODUCTION 



The approximate BCS solution of the pairing Hamiltonian has been extensively used in 



2] . A much 
3j. But, the 



Condensed Matter to study pairing correlations in ultra-small metallic grains [1 
better approximation is given by the Density Matrix Renormalization Group 
pairing Hamiltonian admits an exact solution worked out by Richardson at the beginning 
of the sixties . A more recent derivation of the exact solution can be found in ref. fl. 
The first application of the Richardson exact solution was done in ultra-small grains system 
7, 8|. References 0,0] and jg] marks the resurgence of the Richardson's exact solution of the 
pairing Hamiltonian. The acknowledge to Richardson in refs. j^, s| constitutes a recognition 
to him after forty years in the oblivion. 

The Richardson exact solution has been used to study the effect of the resonant single- 

n n 

particle states on the pairing Hamiltonian [9] . In ref. [10J the authors gave an interpretation 



of the pair energies from the Richardson solution. They relate the pairing correlations with 
the pair energies distribution in the complex plane. The pairing Hamiltonian is not only 



exactly soluble but also quantum integrable llMl4j|. Besides the constant pairing, a very 



special kind of separable pairing interaction also admits an exact solution 



on exact solutions of the pairing Hamiltonian can be found in the ref. 17] 
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Q. 



A review 



The pairing Hamiltonian approximates the influence of the residual interaction acting 
among the valence states lying close to the Fermi level. However, it is an open question 
how one must treat pairing in the continuum. Previous studies on the contribution from the 
continuum to pairing have been reported in refs. Q, [l^ . 

In this paper we reformulate the problem of determining the exact eigenenergies of the 
pairing Hamiltonian when the continuum is included. Real and complex energy representa- 
tions of the continuum are used. The BCS approximation is not a convenient tool to treat 



many-body pairing close to the drip line 20 



2lj . It is the intention of this paper to give 



an exact treatment of the many-body pairing which overcomes the drawbacks of the BCS 
treatment. 

The paper is organized as follows. Section HT1 briefly reviews the derivation of the Richard- 
son equations with the continuum represented on the real energy axis or in the complex 
energy plane. In Sec. |TTT] the low lying states of even Carbon isotopes are evaluated and 
a comparison of the solutions using the real energy representation are compared with the 



ones obtained in the complex energy representation. The trajectory of the pair energies are 
analyzed as a function of the pairing strength. The continuum pair energies are introduced 
in this section. Finally, Sec. [IV] summarizes the main results of the paper. 



II. METHOD 

In this section the Richardson equations for a continuum basis is given. First the contin- 
uum is included by enclosing the system in a large spherical box. After the final equations 
have been obtained, we take the limit of the box to infinity and introduce the single particle 
level density. In order to avoid the Fermi gas we take the derivative of the phase shift for 
the continuum part of the single particle level density 22|. Finally, we parametrized the 
CSPLD for the resonant partial waves and make the analytic continuation to the complex 
energy plane. 

A. System in a Box 

In this sub-section we follow the derivation of the exact solution as it was given by Jan 
Von Delft and Fabian Braun in ref. 6[. The inclusion of the system in a large spherical 
box provides a finite discrete set of negative (bound) energies and an infinite discrete set of 
positive (continuum) energies. Let us called e a the discrete energy with degeneracy 2j a + 1, 
with a = {a,m a } = {n a ,l a , j a ,m a }. The pairing Hamiltonian is given by, 

a am a >0 bmp>0 

with c4 = (— Y a ~ ma c\_ ma . We introduce the pair creation operator 

4 = E c « c * > ( 2 ) 

m a >0 

which creates a pair of time reversal states with quantum number a. 

Following Von Delft and Braun [6|, who were inspired by a suggestion by Richardson, 
we propose the A-body (A = 2A pail .) eigenfunction as the antisymmetrised product of A pair 
wave functions as, 

AT f 

At 

|0> , (3) 




2e n — E n 

t=l \ a 



where the energies E Vi are related to the eigenvalues E of the Hamiltonian Hp by 

N ■ 

1 ' pair 

E=J2 E Pi . (4) 



In order to meet the eigenvalue equation Hp\ty) — E |\&), the parameters called pair 
energies, must verify the following set of A^ pair couple system of equations {f| 

N ■ 

2 „ 2£ « - L Pi 7^ ^ - ^ 

where the first summation contains negative and positive energies. The interpretation of 
this set of equations, called Richardson equations, is that the many-body fermions with 
pairing force behave like the many-boson system with one-body force. Both systems are 
described by the same wave function with the difference that the fermions have to satisfy 
the Richardson equations (|SJ) in order to fulfill the Pauli principle jl, 5]. 



B. Continuum Real Energy 

In making the limit of the box to infinity the single particle states becomes more and 
more dense. In that limit the sum becomes an integral, i. e., 

/oo 
m de . (6) 

The single particle density g(e) is the sum of the bound (negative energy) states plus the 
continuum (positive energy) states. We make the Anzatz that the single particle density in 



the continuum is given by the derivative of the phase shift 



22] 



9(e) = £ (2* + 1) 6{e - e b ) + £ ^ , (7) 

b c 

the index b = (n^, refers to bound states and c = (l c ,j c ) to continuum states. The 
first summation is over the valence bound states while the second one is over the continuum 
partial waves. In practical applications an upper limit / max is set for the number of partial 
waves. 

The Richardson equations in a representation which includes the continuum becomes, 

oo 



E n . — En . 
j-ti Pj Pi 



where the factor db = 2% + 1 — 2Nb takes into account the blocking effect of the Nb unpaired 



states 



5]. The CSPLD becomes, 



9ie) = £ 



2j c + 1 d5 c 
7r de 



(9) 



C. Continuum Complex Energy 

The presence of the single particle resonances appear in the CSPLD, as well as in the 
cross sections, as bumps. They correspond to states in the continuum (positive energy states) 
which are well localized inside the nuclear surface for a time greater than the characteristic 
nuclear time 23(. One can thus split the summation in resonant (r) and non-resonant (nr) 



(background) contributions as, 



9(e) = g n Je) +g Bckg (e) , 
2j r + 1 dS r 



9n CB { £ ) 



E 

r 

E 



7r de ' 

2jnr "I - 1 d5 nr 

7r de 



(10) 
(11) 

(12) 



The sin gle particle density for the resonant states at energies e r and widths T r can be 



written as 



9 Re A £ 



) = E 



2j r + 1 



r r /2 



7T 



e - e r 



+ (r r /2) : 



(13) 



The resonant parameters can be represented by a single complex number e r = e r — 
i T r /2 which corresponds to the eigenvalue of the mean- field Hamiltonian with pure outgoing 
boundary condition 25j. By rotating the integration contour of the resonant part of the 
CSPLD to the negative imaginary axis, and applying the Cauchy theorem, one gets the 
Richardson equations in terms of the complex energy states, 



'-IE 

b 



dh 



G 



2 2eb — E Pk 2 



E 



2jr + 1 G 

2e r E Pk 2 jq 



9 CxBckg 



de°^L-° Hde 

2e - E Pk 2 J 2e- iE„ 



e s-^ 

Pk l-Lh, 



where 



S'cxBckg \£) 



E 



2j r + 1 



r r /2 



7r [e — ie r 



(Y r /2f ■ 



(14) 
(15) 
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In an overstatement (the "density" (^CxBckg CcLIl not be denned outside the integral) one 
could say that the background contribution to the Richardson equation has a real part 
coming from the non-resonant scattering partial wave states g Bckg and a complex contribution 
fi'cxBckg which is a remnant of the complex analytic extension from g Rcs . Because the presence 
of the complex energy Gamow states in the second summation in Eq. (|14p . the complex 
contribution of g CxBckg is necessary to make E = E Pi real. In Eq. (TH1) we have assumed 
that there is not blocking effect due to continuum states. 

For the seniority zero case and neglecting the background, Eq. fTH]) reduces to the 
Richardson equations in the Gamow basis introduced in ref. [9[]. In this case the complex 
pairing energies are not complex conjugate to each other, i.e. E = ^\ E Pi may be complex. 



D. Exact Spectrum 

The solution of the Richardson equations ([8]) with the "boundary condition" , 

lim G ^ +E n = 2e Pi , (16) 

and the blocking effect, determine the ground state and the excited state energies of the 
pairing Hamiltonian. 

The 12 C nucleus has three bound configurations (sec. IHI A lj) . The first (1) and second 



(2) configurations can accommodate a single pair, while the third configuration (3) can 
accommodate three pairs. The configurations (1), (2), and (3) are related to the single 
particle states Opi/2, Is 1/2, and 0^5/2, respectively. Then e Pl = £o Pl/2 , £ P2 = £ 0s 1/2 , and 
e pa = e Pi = e P5 = £od 5/2 - From the bound configurations we can accommodate up to five 
pairs ( 22 C). Because the inclusion of the continuum we will be able to go beyond the nucleus 
22 C. 



1. Ground State 

The ground state (g.s.) configuration for a system with iVp air corresponds to fill the lowest 
-Wpair configurations by solving the Richardson Eq. ([8]) with the blocking coefficient d b = 
2j b + 1 because the g.s. has seniority zero and there are no unpaired states (all N b = 0). For 
example, the g.s. of the isotope 14 C corresponds to solving one single Richardson equation 
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(JHJ) with the boundary condition limc^Q+E Pl = 2e Pl . Let us called this configuration (l) 2 . 
The g.s. of the isotope 16 C corresponds to solving two Richardson equations (jSJ) with the 
boundary conditions lim G ^ +E Pl = 2e pi and lim G ^. +E p2 = 2e P2 . This is the configuration 

(1) 2 (2) 2 , and so on. The ground state energy E is given by Eq. (j3J) with A^ pair = 1,2 and so 
on. 

2. Excited States 

We have to distinguish between excited states with seniority zero and seniority two. 
Seniority Zero (v = 0): The seniority zero excited states are found by solving as many 
equations (jHJ) as pairs, like for the g.s., but with a boundary condition other than the ground 
state. For example, the first and second + excited states of 14 C are found as the solution of 
a single equation with the boundary conditions lim G _> +E p2 = 2e P2 , and lim G ^. + E P3 = 2e P3 , 
respectively. We called such configurations (2) 2 and (3) 2 . As a second example let us consider 
the first + excited state of 18 C. It is found by solving three equations (JSJ) with the boundary 
conditions lim G ^.o+E Pl = 2e Pl , lim G ^ +E P2 = 2e P3 , and lim G ^ +E P3 = 2e P3 . We called this 
configuration (1) 2 (3) 4 . The energy E of the v — excited state is like Eq. (j4]) but using the 
excited pair energies. 

Seniority Two (v = 2): The seniority two states are found by solving iVp air = (A — 12) — v 
equations (jSJ), where A is the mass number of the isotope. This is one equation less than the 
number of pairs in the ground state. The factor db in Eq. (|SJ) is given by db = 2% + 1 — 2Nb, 
where b labels the blocking configuration. For example, to find the v — 2 states in 14 C one 
does not need to solve any equation since iV pa i r = (14 — 12) — 2 = 0. The u = 2 state 
energy is just the sum of the single particle energies E = ei + e m of the unpaired levels I 
and m. Let us assumed that the blocking states for the isotope 16 C are the configurations 

(2) and (3), i. e. N\ = 0, and N 2 = N 3 = 1. Then, we have to solve a single equation with 
d\ — 2, d 2 = and d% = 8 and the boundary condition lim G ^ +E Pl = 2e Pl . Let us call this 
configuration (1) 2 (2)(3) which gives the degenerate levels 2 + , 3 + . The energy of such a state 
is E = E pi +62+ £3. As the last example, let us consider the first 2 + , 4 + states in 18 C. This 
level in found by solving two equations with the boundary conditions lim G _> +E Pl = 2e Pl 
and lim G ^Q+ E p2 = 2e P2 , and with Ni = N 2 = and N 3 = 2. The energy of this last state 
is E = E P1 + E P2 + 2e 3 . 
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E. Determination of the Pairing Strength 



In order to determine the strength G we consider the neutron pairing energy PExp(2iVp a i r ) 



for a system of N = 2iYp air valence neutrons 



26 1 



P Exp (2iV pair ) = 2£(2iV pair - 1) - £(2iV pair ) - £(2iV pair - 2) . 



(17) 



The pairing energy in the Richardson model is related to the last pair energy E PN as 
5], 



follows 



-pRich — 2e Pjv 



Re 



pair ) 



(18) 



By imposing the condition P Exp = Pai C h one finds the strength G which reproduces E. 



Determination of the Resonant Partial Waves 



The criterion to decide whether a given partial wave is resonant is to search for the 
poles Eij = eij — i~y of Sij. A physical resonance should satisfy that the half-life calculated 



with the imaginary part Tu/2 of the pole r 



Kin 2 



is bigger than the characteristic time 



t c = 2.6 x 10 23 x A 1//3 sec. 23[. The physical meaning of this criterion is that the particle 
has enough time to interact with the system before it decays. 



III. APPLICATIONS 
A. Parameters 

This sub-section aims to define the real and complex single particle representations. The 
parameters for the interaction are also set up here. The real energy representation consists 
of a finite discrete set of bound states plus a positive real continuum set of scattering 
states. While the complex energy representation consists of a finite discrete set of bound and 
Gamow states plus a complex continuum set of "scattering states" . In the complex energy 
representation we named resonant continuum the set of Gamow states and non-resonant 
continuum to the scattering states with complex energy. 
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1. Single Particle Representation 



The experimental single particle energies in C were taken from ref. |27|: Eq Pi 



1/2 



-4.946 



MeV, e ls = -1.857 MeV, and e od = -1.093 MeV. The single particle density of 13 C was 



calculated with the program 



28 



MeV, V so = 10.5 MeV, a 



with the following Woods-Saxon parameters: V = 55.1 
0.7 fm, tq = r so = 1.27 fm. Fig. [1] compares the CSPLD 



for Z max 
system. 



10 and L 



15. It shows that a cut-off of Z = 10 in Eq. (Q is enough for this 




10 15 20 
E(MeV) 



30 



FIG. 1. Neutron CSPLD in 12 C for two different angular momentum cutoff. 



The negative contribution in Fig. [T]is due to the dominance of the S\/2 state at low energy. 
In accordance to the Levinson theorem, it must be a negative contribution for each bound 
state. For the S\/2 state this negative contribution is close to the continuum threshold. The 
resonant behavior around 2 MeV is due to the resonant state dy?., while the one around 
10 MeV is due the wide resonance fa/2. Using the code Gamow 28] we find the following 
energies for these two states, e 0d , j/2 = (2.2671; -0.416) MeV, and e / 7/2 = (9.288; -3.040) 
MeV. 



2. Pairing Strength 

From the experimental mass excess table we got for the pairing energy Pexp of the isotope 
14 C, -Pexp( 14 C) = 1.516 MeV. In the Richardson model the pairing energy is related to the 
pair energy through P(2V pair ) = 2e PN . — E PN (Sec. Ill Ej) . For 14 C, iVp air = 1, then 



* pair c -"pair 



P( 14 C) = 2e Pl - E P1 with 2e Pl = -9.989 MeV and E Pl = -11.408 MeV. In order to 
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reproduce E Pl with a cutoff energy at 30 MeV, one must take G = 0.7786 MeV. Using the 



parametrization G = ^ we obtained x = 10.900 for A = 14. This value of % is used for all 
Carbon isotopes. Table fl] lists the value of the pairing strength for each Carbon isotope. 

TABLE I. Pairing strength used for the Carbon isotopes. 

Isotope G[MeV] 

14 C 0.7786 

16 C 0.6813 

18 C 0.6056 

20 C 0.5450 

22 C 0.4955 

24 C 0.4542 



B. Results: Real Energy Representation 

After the model space and the interaction are set up one can evaluate physical magnitudes. 
In this subsection we are going to calculate the ground state energy of the carbon isotopes 
14 C to 24 C and the low energy spectrum of the isotopes 14 C to 20 C. 



1. Ground-state Energy 

Solving the Richardson equations ([H]) for the ground state of each carbon isotope, we 
obtained a set of pair energies E± (we set E { for E p .) as it is shown in table [III Complex 
pair energies appear in complex conjugate pairs to give a real eigenenergy. The distribution 
of the pair energies gives information about the structure of the many-body wave function. 
As the many-body state becomes more collective, more pairs accommodate themselves in a 



parabola-like distribution (lOj. Let us quantized roughly the degree of collectivity 7 as the 
ratio of the number of pairs which participate in a parabola versus the total number of pairs. 
We will do this for system with at least four pairs. We observe a high degree of collectivity 
as one approaches the threshold, while the collectivity abruptly drops in the continuum. 
Figures El and H] show the distribution of the pair energies in the complex energy plane 
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for the isotopes 20 C, 22 C and 24 C, respectively. 



TABLE II. Pair energies Ei and ground state energies Eq relative to carbon C for the Carbon 
isotopes 14 C— 24 C. We used Ei for E Pi . The collectivity parameter 7 was defined in the text. 



Isotope 


A'pair 


Ei[MeV] 


E [MeV] 


7 


14 C 


1 


E 1 = -11.398 


-11.398 




16 C 


2 


E x = -10.681 
E 2 = -6.370 


-17.051 




18 C 


3 


E x = -10.495 

E 2>3 = (-4.950; ±1.262) 


-20.394 




20 C 


4 


Ex = -10.379 

E 2 = -4.502 

E 3)4 = (-3.667; ±1.546) 


-22.194 


0.75 


22 C 


5 


Ex = -10.302 

E 2>3 = (-3.729; ±0.110) 

£ 4 5 = (-2.578; ±1.361) 


-22.915 


0.8 


24 C 


6 


Ex = -10.254 

E 2 = -3.924 

E 3 = -3.099 

E 4>5 = (-2.479; ±0.969) 

E 6 = 2.630 


-19.605 


0.5 



Table HIl also shows the ground state energy E of the Carbon isotopes 14 C to 24 C. Fig. [5] 



compares the calculated ground-state energy with the experimental one 29] . It is found that 
the exact solutions follow the overall trend, i.e. the binding energy decreases faster at the 
beginning of the chain and decelerates when it approaches the drip line. The agreement with 
data worsen as the number of neutrons increases. Even when the pairing interaction is a 
schematic one, and not realistic, this investigation suggests that the nucleus 24 C is unbound. 
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FIG. 2. Distribution of the four pair energies in 20 C isotope (dark dots). The white diamond 
correspond to the pair energies for G = 0, i.e. Ej, = 2ej. 
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FIG. 3. Like fig. [2] for the five pair energies in 22 C. 



2. Carbon Isotopes Spectrum 

It is worthwhile to compare the experimental spectrum with the exact solutions of the 
schematic pairing Hamiltonian corresponding to the cases of seniority-zero and seniority-two. 

14 C Spectrum: Table III II gives the excitation spectrum (last column) with respect to 
the ground state configuration (l) 2 . The seniority u, the pair energies and the number of 
pair iVp air = (A — 12) — u (A the mass number) are also given. Figure E] compares the 
calculated levels from table IIHI with that of the experimental one. The quantum number 
of the first excited state 1~ is correctly found with 1.5 MeV less energy. The 0^" and 0% 
excited states are underestimated with respect to the experimental one by 0.951 MeV and 
1.517 MeV respectively. The state 3~ is found 1.37 MeV below the experimental one. The 
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FIG. 4. Like fig. [2] for the six pair energies in 24 C. 
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FIG. 5. Carbon isotopes ground-state energy. 

splitting between the states 3~ and 0% is well reproduced: 280 keV versus the experimental 
175 keV but in inverse order. We missed the first 2 + state and found a 2 + at only 129 
keV from the second experimental 2 + . The near degenerate experimental states 2 + and 4 + 
around 10 MeV are well reproduced. 

16 C Spectrum: Table IIVI shows the pair energies and the excitation spectrum with 
respect to the ground state configuration (1) 2 (2) 2 . Fig. [7] compares the calculated versus 
the experimental spectrum of 16 C. The first excited 2 + state does not appear in our spectrum. 
The first + excited state is very well reproduce with a difference of only 21 keV. We found a 
2 + state at 3.274 MeV which may correspond to the experimental 2 state at 3.986 MeV. The 
first 4 + excited state is found only 125 keV below the experimental one. The experimental 
(3~) is 406 keV from the 3~ calculated state. In the exact spectrum appears a third + state 
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TABLE III. Excited and 


pair enerj 


$ies of 14 C. The enerj 


pes are in MeV. 




Config 






N- pair 






Ex 


(I) 2 


n 
u 


u 


1 


£ pi = -11.398 


-11.398 





(1)(2) 


o 

z 


u , 1 







-6.803 


4.594 


(1)(3) 


2 


2~,3- 







-6.039 


5.358 







0+ 


1 
1 


-C/p 2 — — D. / DU 


— 0. ( DU 


D.Doo 


(3) 2 





0+ 


1 


£ P3 = -3.168 


-3.168 


8.229 


(2) (3) 


2 


2+ 







-2.950 


8.447 


(3) (3) 


2 


2+,4+ 







-1.093 


10.304 



10304 2 , 4 
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8229" 
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FIG. 6. Exact low energy spectrum of C for seniority zero and two compared with experimental 



levels 



The energies are in keV. 



which does not appear in the experimental spectrum. Finally, the expaerimental (4 + ) state 
is 938 keV from the 4 + calculated one. Summing up what we found for the nucleus 16 C, the 
first + , 4 + and 3~ are reasonable well described by the pairing interaction. 

18 C and 20 C Spectra: Tables [V] and [VI] show the pair energies and the excitation 
spectrum with respect to the ground state configuration for the three and four pair systems 
18 C and 20 C respectively. Figure [S] shows the calculated exact eigenvalue of the pairing 
Hamiltonian for 18 C and 20 C. Experimentally, only one excited state in 18 C is known. It 
is a (2 + ) state at 1620 keV from the (0 + ) ground state. Considering what we learn in the 
previous spectra one may place some confidence on the theoretical estimation for the levels 
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TABLE IV. Like table [III] for 16 C. 



Config 


V 


State 


jV air 


EfefMeVI 


ElMeV] 

L J 


ExlMeV] 

L J 


(1) 2 (2) 2 





0+ 


2 


E 1 = -10.681 


-17.051 













E 2 = -6.370 






(1) 2 (3) 2 





0+ 


2 


Ex = -10.823 


-14.003 


3.048 










E 3 = -3.180 






(1) 2 (2)(3) 


2 


2+,3+ 


1 


E 1 = -10.827 


-13.777 


3.274 


(1) 2 (3)(3) 


2 


2+,4+ 


1 


E 1 = -10.848 


-13.034 


4.017 


(2) 2 (1)(3) 


2 


2-,3~ 


1 


E 2 = -5.309 


-11.348 


5.703 


(2) 2 (3) 2 





0+ 


2 


E 2>3 = (-4.777, ±1.079) 


-9.554 


7.497 


(2) 2 (3)(3) 


2 


2+,4+ 


1 


E 2 = -4.827 


-7.013 


10.038 



10038 



7497 



5703 



2,4 



2, 3" 



.(2) 
(3)= 

5^ 



(2 , 3') 



10390 
) 

9420 



19980 
i 

1-9100 
N 8920 



7740 



6109 
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3274 



3048 
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3027 



4088 
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1766 



0.0 



0.0 



FIG. 7. Like fig. [6] for 16 C. The experimental levels above 6.11 MeV are from 
are in keV. 



30l | . The energies 



0+, 4+ and 3f . 



C. Results: Complex Energy Representation 



The first step in the determination of the complex representation is to find the resonant 
partial waves. This is done by evaluating the outgoing solutions (Gamow states) of the 



Schrodinger equation [25|, 



311 ] of the mean field Hamiltonian defined in Sec. IIII A II Then, 
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TABLE V. Like table HI] for 18 C. 



Config 


V 


State 


j\T air 




' L J 


ElMeV] 

L J 


ExlMeV] 

L J 


(1) 2 (2) 2 (3) 2 





0+ 


3 


-^2 3 


£i = -10.495 

= (-4.950; ±1.262) 


-20.394 





(1) 2 (2) 2 (3)(3) 


2 


2+,4+ 


2 




Ex = -10.531 
^2 = -4.809 


-17.525 


2.869 


(1) 2 (3) 2 (2)(3) 


2 


2+,3+ 


2 




Ex = -10.543 
£ 3 = -3.710 


-17.203 


3.191 


(1) 2 (3) 4 





0+ 


3 


-^2,3 


Ex = -10.549 

= (-3.236; ±0.474) 


-17.020 


3.374 


(1) 2 (3) 2 (3)(3) 


2 


2+,4+ 


2 




Ex = -10.573 
E 3 = -2.630 


-15.388 


5.006 


(2) 2 (3) 2 (1)(3) 


2 


2 _ ,3- 


2 


-^2,3 


= (-4.177; ±0.772) 


-14.393 


6.001 


(3) 4 (1)(2) 


2 


o-,i- 


2 


-^2,3 


= (-3.576; ±0.981) 


-13.955 


6.439 


(2) 2 (3) 4 





0+ 


3 


-^1,3 


E 2 = -4.405 
= (-3.627; ±1.433) 


-11.659 


8.735 


(2) 2 (3) 2 (3)(3) 


2 


2+,4+ 


2 




E 2 = -4.008 
E 3 = -3.243 


-9.347 


11.047 



11047 



2\4 T 



6001 



3191 



8735 


+ 


6439 


0", f 


9 


5006 


2 + ,4 + 


3374 


+ . 




2869 


2 + ,4 + 


0.0 


+ 



0, 1 



6993 



2. 3 


6225 


2 + ,3 + 


3147 


2 + .4 + 


2381 


+ 


0.0 



18, 



20, 



FIG. 8. Exact low energy spectra of 18 C and 20 C for seniority zero and two. The energies are in 
keV. 
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TABLE VI. Like table QH for 20 C. 


Config 


;/ 


State 


7V Dair EAMeV] 


E[MeV] 


Ex[MeV] 


(1) 2 (2) 2 (3) 4 
\ ± J \ u ) 





0+ 


4 Ei = -10 379 


-22.194 











E 2 = -4.502 












E 3A = (-3.667; ±1.546) 






(1) 2 (2) 2 (3) 2 (3)(3) 


2 


2+,4+ 


3 Ei = -10.369 


-19.813 


2.381 








E 2 = -3.993 












E 3 = -3.238 






m 2 (3) 4 (2H3) 


2 


2+ 3+ 


3 Ei = -10.406 


-19.047 


3.147 








E 3>4 = (-2.846; ±0.665) 






(2) 2 (3) 4 (1)(3) 


2 


2-3- 


3 E 2 = -4.080 


-15.969 


6.225 








E 3A = (-2.925; ±0.896) 






(3) 6 (1)(2) 


2 


o-,i- 


3 E A = -3.167 


-15.201 


6.993 








E 2t3 = (-2.615; ±1.113) 






the half-life of the Gamow state is compared with the characteristic time of the system 


t c = 5.953 x 1(T 23 


sec 


fsee Sec. IIIFD. Table IVIII compares the characteristic time with the 


half-life of the states £od 3/2 an d 


£o/ 7 / a - The half-life of the state Od 


3/ 2 is around nine times 


bigger than the characteristic time. The Of 7/2 state seems to be a 


wide resonance, but the 


comparison with the characteristic time shows that its half-life is a 


bit bigger 


than r c . 


TABLE VII. Comparison of the half-life versus the characteristic time (Sec. 


IIIFD. 


state 






T 1/2 [sec] 




Ti/ 2 /t c 








5.485 x 10" 22 




9.21 


0/ 7 /2 






7.505 x 10~ 23 




1.26 



The effect of the resonant continuum was already investigated in ref. {q|. In order to 
investigate the effect of the non resonant continuum on the many-body correlations we 
compare in fig. [9] the ground state energy of the nucleus 22 C as a function of the pairing 
strength for three different model spaces: (i) Bound: {0pi/ 2 ,lsi/ 2 ,0d 5 / 2 }, (ii) Resonant: 
{Opi/2, lsi/2, Ods/2, 0d 3 / 2 , O/7/2}, an d (in) Continuum (Secc. IIII Bj) . It is observed that the 
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resonant and non resonant continuum states can be neglected as long as the pairing force 
is not very strong 9|. As the interaction increases the continuum starts to be important. 
The curve labeled as "Continuum Representation" gives the ground state energy when the 
resonant and non resonant continuum is included in the representation through the CSPLD. 
The figure shows clearly the energy gain due to the inclusion of the continuum. The curve 
labeled as "Resonant Representation" gives the energy when only the resonant states are 
included in the representation. For very big strength the non resonant continuum becomes 
as important as the resonant continuum. 




Let us compare the evolution of the pair energies Ei in the bound and the resonant 
representation versus the pairing strength. Figure [TU1 shows Ei for G from G = 1.0 MeV 
to G = 0.005 MeV in the nucleus 22 C . The continuum (dot) line corresponds to bound 
(resonant) representation. The deeper pair energy E\ is little affected by the model space 
(one can not distinguish between the two curves). The other pairs are more affected for 
big value of the strength. The difference diminishes as the interaction decreases. The same 
effect was observed in the ground state energy (fig. [9]). The pairs E 2 and E 3 are complex 
conjugate partners for G > 0.51 MeV and they move at the same pace as G changes. When 
they become real E 2 approaches to the uncorrelated pair energy 1e 2 while E% moves faster 
to the uncorrelated pair energy 2^3. The pairs E4 and E$ remain complex conjugate for all 
no zero values of the strength. 

As a last application we will calculate the evolution of pair energies in the continuum, 
i.e. pair energies with positive real component. To this aim let us study the nucleus 28 C 
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-10 -9 -8 -7 -6 -5 -4 -3 -2 
Re[E.] (MeV) 



FIG. 10. Pair energies in the ground state C versus pairing strength G for G = 1.0 MeV to 
G = 0.005 MeV. The continuum line corresponds to the bound representation while the dot line 
corresponds to the resonant representation. The arrows point in the direction of decreasing G. 

with eight pairs. Fig. [TT] shows the evolution of the pairs for strength from G = 2.2 MeV 
to G = 0.2 MeV. The bound (negative real component) pairs E\ to E 5 follow a trajectory 
similar to that in 22 C with the difference that the complex partners E2 — E% and E4 — E§ are 
only approximately complex conjugate to each other. They become truly complex conjugate 
partners as the interaction approaches zero. On the other hand, the pairs in the continuum 
show a striking behavior. The typical movement to the right is not follow by all the positive 
energy pairs, i.e the continuum pairs may converge to its uncorrelated energy from right or 
left as G decreases. Besides, the pairs seem to converge to the real part of the uncorrelated 
pair energy lim G ^ Q +Ei = 2Re[ei\ when £j is a Gamow state. 



IV. CONCLUSION 

The contribution of this paper to the exact solution of the pairing Hamiltonian is the 
inclusion of the resonant and non resonant continuum through the continuum single particle 
level density (CSPLD). The Gamow states, which appear in the complex energy representa- 
tion, provide the main contribution from the continuum. It is worthwhile to point out that in 
the representation these states have exactly the same status as bound states. The difference 
is that the states in the continuum are no affected by blocking effects. The inclusion of the 
continuum has allowed us to study the unbound isotope 24 C and beyond. It was found that 
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FIG. 11. Evolution of the pair energies in the ground state of C as a function of the pairing 
strength G from G = 2.2 MeV to G = 0.2 MeV. The arrows point in the direction of decreasing G. 

the continuum pairs (pair energies with positive real components) converge to the real part 
of the uncorrelated pair energy and they do not appear in complex conjugate partners. As 
a consequence the total energy may be complex. It was shown that from the exact solution 
of the pairing Hamiltonian the CSPLD can be used to investigate the effects of the resonant 
and non resonant continuum states upon the many-body pairing correlations. 
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